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OVERVIEW: UPPER STAGE MODELING 
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http://www.boeing.com/defense-space/space/delta/delta4/d4h_demo/book04.htrnl 
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OVERVIEW: WHAT CAN HAPPEN INSIDE TANKS? 
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t = 1 ,229.5 sec 




Parking orbit 
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Transfer orbit 
5.14-hr coast 


SECO-3 

t = 20,427.3 sec 

DemoSat payload separation 

t = 20,977.5 sec 

Orbit = 19,622 x 19,623 nmi 

at 10-deg Inc 


Stage exposed to solar heating 

Propellants (LH 2 and LOX) may thermally 
stratify 

Propellants may boil 

Slosh events during maneuvers 

Upper stage must re-start at conclusion of 
coast phase for insertion 



http://www.boeing.com/defense-space/space/delta/delta4/d4h_demo/bookl4.html 
XSS-10 view ,ot Delta II rocket: An Air Force Research Laboratory XSS-10 micro-satellite uses its onboard camera system to view the second stage of the 
; Delta II rocket during mission operations Jan. 30. (Photo courtesy of Boeing.), http://www.globalsecurity.org/space/systems/xss.htm 
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Engine Inlet Temperature 


OVERVIEW: WHAT CAN HAPPEN INSIDE TANKS? 


• Propellant T&P must be within specified range for turbomachinery operation 

- If propellants outside specified T&P box engine may not restart 

- Orbit cannot be circularized 




.pratt-whitney.com/prod_space_rl 1 0.asp 
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MOTIVATION 


Rotation present during missions to evenly heat spacecraft 

Effect rotation has on propellant thermal properties unknown 

Upgrade current analytical/numerical stratification models to include rotation 
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MISSION PARAMETER RANGES 

• Tank Dimensions: 

- Square 3 m diameter tanks 

• Cryogenics: LH 2 , LOX 

- T bulk LH 2 : 16 K, 28.8 °R, -430.9 °F 

- T bulk LOX: 91 K, 163.8 °R, -295.9 °F 

• Tank Pressure (All Cryogenics): 30 psi 

• Initial Fill Levels: 10, 20, 30% 

• Heating Conditions: 

- Constant wall temperature: 0 = T wall - T bulk : AT = 0. 1 , 0.5, 1 .0 K 

- Heat flux to fluid: 5-100 W/m 2 

• Reduced Gravity Environment: g/g 0 = 10" 4 , 10' 3 , 10' 2 , 10 1 , 1 

• Rotation rates: co = 0.1, 1,5 °/sec 

• Orbital Transfer Time (Simulation Time): 2-4 HR 
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MASTER MODEL: BASIC FRONT END OPTIONS 


1. Tank geometry 

- Tank diameter, height 

- Square bottom 

2. Boundary layer nature and heat transfer coefficient selection 

- Free convection 

• Laminar/Turbulent (w/ & w/out switching) 

3. Wall temperature settings 

- Constant inner wall temperature 

- Constant inner wall heat flux 

4. Rotation rate 

5. Gravity level 
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GENERAL MODELING PHILOSOPHY 
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• Stratum growth A(t) 

- u(y) depends on if heating is constant wall 
temperature or constant heat flux, q 

- u(y) depends on nature of boundary layer 

- Provides differential equation for A(t) 

• Stratum temperature, T s (t) 

- Heat entering side wall into boundary layer 
is used to increase stratum temperature 

- Energy exchange with ullage negligible 

- T s assumed uniform 
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RELEVANT NON-DIMENSIONAL NUMBERS 

• Grashof number, Gr, governs heat transfer regime for 
3 constant wall temperature 

gpOL 

- Ratio of buoyancy to viscous forces 

- p, Volumetric thermal expansion coefficient 

- 0, Wall to Bulk temperature difference 


V 


Ra = Gr Pr 


• Rayleigh number, Ra, is product of Grashof and usual 
Prandtl number, Pr 

• Prediction of boundary layer transition 

- If Ra < 1 0 9 — > Laminar 

- If Ra > 10 9 — > Turbulent 



Modified Grashoff number, Gr*, governs heat transfer 
regime for uniform heat flux, q w 


Ra = Gr * Pr 


Modified Rayleigh number, Ra*, for uniform heat flux 


• Others: Reynolds Number, Re (momentum to viscous) Weber number, We (inertial to 
capillary), Froude number, Fr (inertial to body), and Bond number Bo (body to capillary) 
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Ra and Ra* vs. g/g 0 MAPS for LH 2 and LOX 


R vs. g/g (LH2) 
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Maps laminar or turbulent boundary layers possible for typical mission profiles (NIST data) 
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WHAT IS IMPACT OF BBQ ROLL ROTATION: 

Typical rotation rate, co~l°/sec Shape at g/g 0 =l, co~850°/sec 



• Assume liquid is in solid body rotation (transients can also be treated) 

• Model extra height that liquid gains along wall as a longer interfacial heat transfer length 

• Center point in radial direction of tank is taken to be point where percent of bulk remaining 


is referenced — » worst case scenario 

• Trade off between heated area and surface area to distribute warm stratum 
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ROTATION / STRATIFICATION COMBINED MODEL 
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ROTATIONAL CASES 


Re-examine boundary layer / stratum mass balance 
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Re-derive energy balance to take into account additional heating area 
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COMBINED ROTATION / STRATIFICATION MODEL: LH 2 and LOX 


Fraction Remaining vs. Time (LH2, g/g o = 10' 4 ) 



Fraction Remaining vs. Time (LOX, g/g o = 10' 4 ) 



T strat VS - Tlme (LH2> g/g o = 1 0_4) T strat VS Time ( LOX ’ 9^ = 1 °' 4) 



For q=10 W/m 2 , L=3, R=1.5, 20% fill level, H/R=0.4 and co=l°/sec at g/g 0 =10- 4 : 

- Rotation decreases time to stratification time by ~ 15% 

- Rotation increases stratification temperature by ~ 1 .0 K 
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EFFECTS OF ROTATION AND TRADEOFFS 

Increased boundary layer running length (H — > HJ 

- more heated area 

- larger Grashof number 

Larger surface area at bulk-stratum interface (S — > S paraboloid ) 

- increased mass flow rate into stratum layer 

- more area to spread mass flow 




0S (K) 


EFFECTS OF ROTATION 

• Spinning always increases stratification 

• Stratum temperature affected by spin rate; especially at low gravity levels 

• LOX cases shown with heat flux of 5 W/m 2 after 2 hour mission 


Effect of Rotation on Stratum Temperature 
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CO critical (deg/s) 


EFFECTS OF ROTATION 


• ® critical spin rate to minimize stratum temperature 

• co critical needed for large g/g 0 impractical 

• (^critical ne ^ded for typical mission profiles very practical (co < 1.5 deg/s) 

• LOX results discussed previously shown 



— g_go = 1e-3 
— g_go = 1e-4 





SUMMARY/ CONCLUDING REMARKS 

Thermal stratification impacts T&P at conclusion of coast phase 

Rotation (creeping of fluid up side walls) has large effect for co=17s and g/g 0 =10 -4 

- ‘Classical’ literature model upgraded to include rotation effects 

- Can decrease time to stratify by 30-60 minutes during 4 hour coast 

- Larger heating area and lower liquid level above sump inlet 

- For various missions stratum temperature may increases or decrease relative to 
case 

- Mixed tank temperatures always larger because A increased with rotation 
Future work 

- Comparison with CFD studies 
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